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Abstract 



^ \ We show that the gravitino-overproduction problem is prevalent among inflation 

models in supergravity. An inflaton field cj) generically acquires (effective) non- 
vanishing auxiliary field Q^^^^ , if the Kahler potential is non-minimal. The infiaton 
field then decays into a pair of the gravitinos. We extensively study the cosmological 
constraints on Q^^^ for a wide range of the gravitino mass. For many inflation 

models we explicitly estimate Gi^^\ and show that the gravitino-overproduction 
problem severely constrains the inflation models, unless such an interaction as A' = 
k/2 \4>\'^z^ -|-h.c. is suppressed (here z is the field responsible for the supersymmetry 
breaking). We find that many of them are already excluded or on the verge of, if 
K~ 0(1). 



1 Introduction 



The gravitino is the most important prediction of unified theory of quantum mechanics 
and general relativity such as the superstring theory (i.e. supergravity (SUGRA) at 
low energies) |1 . However, the presence of the gravitino leads to serious cosmological 
problems depending on its mass and nature. If the gravitino is unstable and has a mass 
m3/2 in the range from O(IOO) GeV to 0(10) TeV, the decay of the gravitino destroys 
light elements produced by the big-bang nucleosynthesis (BBN). To keep the success of 
BBN the reheating temperature Tr after infiation should be lower than O{10^~^) GeV 
suppressing the gravitino production by thermal scattering. On the other hand, if the 
gravitino is light as 1713/2 < 0(10) GeV and it is stable (that is, the lightest supersymmetric 
particle (LSP)), the reheating temperature should satisfy Tr^ 0(10'') GeV(m3/2/l GeV) 
for m3/2 ^ 100 keV for the gravitino density not to exceed the observed dark matter 
density. 

In a recent article j2], we have pointed out that there is a new gravitino problem 
beside due to the thermal production of the gravitino. That is, an infiaton field has 
nonvanishing supersymmetry(SUSY)-breaking auxiliary field (or more precisely Q^^^ 
as will be defined later) in most of infiation models in SUGRA, which gives rise to an 
enhanced decay of the infiaton into a pair of gravitinos, if the Kahler potential is non- 
minimal. Thus, we have stringent constraints on the (effective) auxiliary field Q^^^ to 
suppress the production of gravitinos in the infiaton decay [2^ . This gravitino production 
in infiaton decay is more effective for lower reheating temperature, while the production 
by particle scatterings in the thermal bath is more important for higher temperature Tr. 
Therefore, the direct gravitino production discussed in this paper is complementary to 
the thermal gravitino production, and the former may put severe constraints on infiation 
models together with the latter. 

The purpose of this paper is to discuss this new gravitino problem in a broad mass 
range of the gravitino including 7713/2 — 0( 100) TeV region suggested from anomaly- 
mediated SUSY breaking models jS]. We assume, in the present analysis, that there is no 
entropy production after the end of reheating by the infiaton decay. However, we briefiy 
discuss, in the last section of this paper, the case that a late-time entropy production 
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takes place. 

In Sec. 121 we briefly review the gravitino problem in cosmology and in Sec. El we 
calculate the abundance of gravitinos produced by particle scatterings in the thermal 
bath and show cosmological constraints on the reheating temperature Tr. In Sec. |31 we 
discuss the enhanced decay of the inflaton into a pair of gravitinos and give cosmological 
constraints on the (effective) auxiliary fleld G^^\ In Sec. El we explicitly calculate the 
precise value of Q^^^ for inflation models in SUGRA to demonstrate how severe the new 
constraints are. The last section is devoted to conclusions. 

2 Gravitino problem 

The gravitino is the SUSY partner of the graviton in SUGRA and it acquires a mass 
in a range of (9(100) GeV — (9(10) TeV in gravity-mediated SUSY-breaking models 
Such a gravitino is likely unstable and its lifetime is very long because interactions of the 
gravitino are suppressed by inverse powers of the reduced Planck scale Mp. The grav- 
itino dominantly decays into the standard-model (SM) particles and their superpartners, 
which may produce a large entropy and destroy the light elements synthesized in BBN. 
As a result, the predictions of BBN may be signiflcantly changed unless the primordial 
abundance of the gravitino is sufficiently small 

In gauge-mediated SUSY-breaking models [7^, the gravitino is light (m3/2 ^ 10 GeV) 
and stable. In this case the gravitino may give too much contribution to the present 
cosmic density of the universe. 

In the inflationary universe, the primordial gravitino is once diluted but it is produced 

during reheating epoch after the inflation. Thus, even in the inflationary models, we may 

still have the gravitino problem [Hj. As shown in the next section, this leads to very 

stringent constraints on the reheating temperature Tr since the gravitino abundance is 

approximately proportional to Tr. The constraints are given in Pl inH ITT IIT^ITHllT^ITHl lTBj 

for the unstable gravitino and in J7] for the stable one. 

#^Although the gravitino mass can be either much hghter 0] or much heavier in no-scale models, 
we do not consider such possibilities in this paper. 
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3 Thermal production of gravitinos and cosmological 
constraints on the reheating temperature Tr 



In this section we show the abundance of the gravitinos thermally produced after inflation 
and derive constraints on the reheating temperature. 

During reheating the gravitino is produced through scatterings of particles in ther- 
mal bath. The interactions of the gravitino with a gauge multiplet (A^, A) and a chiral 
multiplet {rj, x) are described by 

1 „ .V ., .. 1 



V2Mi 



V2Mi 



8M, 



-V^M [7^ 7l t'^AF, 



(1) 



where Fj^p is the field strength of the gauge field. (Here, D^, denotes the covariant derivative 
and xr satisfies (1 — 75)X-R = 0.) The thermally averaged cross section of the gravitino 
production for an SU (N) super Yang-Mills model with rif pairs of fundamental and anti- 
fundamental chiral superfields is calculated in Ref. fH] as 



(cTfrel) 



3/77-3^2 



X 



C(3)^ 



327rM| 

ln(TVm^ + 0.3224] (A^ + m) + 0.5781nf } 



(2) 



where rrig is the gaugino mass and mg,th is the thermal mass of the gauge boson which is 
given as m^^-j^ = {l/6)g'^{N + n{)T'^. 

Solving the Boltzmann equation with the above cross section, one can obtain the 
gravitino-to-entropy ratio Y3/2 which is well approximated by ^IS] 

,2 



3/2 



X 



1.9 X 10"^2 
1 + 0.045 In 




R 



1010 GeV 
- 0.028 In 



1010 GeV 



(3) 
Tr. 



,1010 GeV. 

where we have taken = 3 for QCD and rrig^ is the gluino mass evaluated at T 
Notice that the gravitino abundance is roughly proportional to Tr. 

For the gravitino of a relatively large mass ^ 100 GeV, it likely decays to the SM 
particles and their superpartners. In that case, high energy photons and hadrons emitted 
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in the gravitino decay may destroy the hght elements (D,^He, ^He, ^Li, ■ • ■) and hence 
spoil the success of BBN. Since the gravitino abundance is approximately proportional to 
Tr, we obtain an upper bound on the reheating temperature after inflation. 

Energetic photons from the radiative decay of gravitino {ipf^ — > 7 + 7) deconstruct D, 
which gives an upper bound Tr^ 10^ GeV for 7713/2 ~ 1 — 3 TeV. They also cause an 
overcreation of '^He due to photo-dissociation of ^He, which leads to the most stringent 
constraint on Tr as ^ 10^ - 10^ GeV for mg/a ^ 100 GeV -1 TeV. 

However, it was found in Ref. [19j that the hadronic decay gives a more stringent 
constraint on the abundance of gravitinos and equivalently on the reheating temperature 
Tr because mesons and nucleons produced in the decay and subsequent hadronization 
processes significantly affect BBN. In particular, when the branching ratio into hadrons is 
~ 1 as expected for the gravitino decaying into gluino and gluon, the effect of the hadronic 
decay is much more serious than the radiative one. 

In the case of Bh = 1 {Bh- the branching ratio of the hadronic decay), the upper bound 
on Tr for relatively light gravitino m3/2 ~ 0.1 — 0.2 TeV comes from the overproduction 
of ^He as 

Tr < (1 - 4) X 10^ GeV for 7713/2 - 0.1 - 0.2 TeV {Bh ~ 1), (4) 

Here we conservatively assume rrig -C 7713/2. For 7713/2 — 0.2 — 1 TeV, non-thermal pro- 
duction of ^Li sets the very stringent constraint as, 

Tr < 3 X 10^ -4 X 10*^ GeV for 7773/2 ~ 0.2 - 2 TeV (Bh^l). (5) 

For larger gravitino mass the destruction of D gives the stringent constraint, 

Tjj < 5 X 10^ - 1 X 10^ GeV for 7773/2 -2-10 TeV {Bh ~ 1). (6) 

Since the gravitino of mass larger than 10 TeV decays before the light elements are 
synthesized, the stringent constraint is not obtained from hadro-dissociation processes. 
However, the mesons (mainly pions) produced at ~ 1 sec alter the proton-neutron ratio 
and increase the abundance of ^He, from which the upper bound on Tr is obtained as 

Tr < (3 - 10) X 10^ GeV for 7773/2 -10-30 TeV. {Bh ~ 1) (7) 
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For m3/2 ~ 30 TeV the gravitino decay little affects BBN in the case of -B/i = 1. 

When the main decay mode is not hadronic, the above constraints become milder. 
However, even if the gravitino dominantly decays into a photon and a photino, the 
hadronic branching ratio is non-vanishing since the quark-anti-quark pair can be at- 
tached at the end of the virtual photon line. In this case, Bh is expected to be of order 
C^(aem/4vr) ~ 10^3. Even such small B^ makes the constraint severer than that for pure 
radiative decay {B^ = 0) as 



Tr < 1 X 10*^ - 3 X 10^ GeV for 
Tr < 1x10^ -3x 10^ GeV for 
Tr < 2 X 10^ - 1 X 10^ GeV for 



7113/2 -0.1-1 TeV (8) 

m3/2 - 1 - 3 TeV (9) 

m3/2 - 3 - 10 TeV (10) 
(Bh ~ 10-=^) 



where the upper limits on Tr are imposed by ^He overproduction, ^Li overproduction and 
D destruction, respectively. In the case of B^ = 10~^, no sensible BBN bound exists for 
m3/2 > 10 TeV. 

The corresponding constraints on Y3/2 which will be used later are obtained by sub- 
stituting the upper bounds on Tr into 0, 
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30(10) TeV with Bh 


= 1(10-3) 
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(12) 



straints are obtained from BBN. However, another constraint comes from the abundance 
of the LSP produced by the gravitino decay. Since the gravitino decay temperature is 
rather low, one LSP remains as a result of the decay of one gravitino. The relic LSP 
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density is 

where mLSP is the LSP mass, and we have conservatively neglected the contribution from 
the thermally produced LSPs. According to the recent WMAP result 20j, the dark 
matter density is Qj^uh"^ — 0.11 ± 0.01 {h: Hubble parameter in units of lOOkm/s/Mpc). 
Requiring the LSP density smaller than the upper bound on the dark matter density at 
95 % C.L., we obtain 

which is applicable for the unstable gravitinos. This bound is important especially for 
the gravitino heavier than 30(10) TeV, which falls in the range suggested from anomaly- 
mediated models of SUSY breaking ^3^. In the anomaly-mediated SUSY breaking models, 
the LSP is mostly composed of the wino W and its mass is related to the gravitino mass 
as 

= — m3/2 ~ 2.7 X 10"^m3/2, (15) 

92 

where P2 and (72 are the beta function and the gauge coupling of SU{2)l. Since the 
thermal relic of the wino LSP is less than the observed dark matter abundance as long as 
my^r^2 TeV EB, we obtain 

Tr < 9.3 X 10^ ( '^%\, ] ' GeV, (16) 
" ~ VlOOTeV/ ^ ^ 

for m3/2 ^ 7 X 10^ TeV in the anomaly-mediated SUSY breaking models with the wino 

LSP. 

When the gravitino is light (^ 10 GeV) which is expected in gauge-mediated SUSY 
breaking models, the gravitino may be the LSP and hence stable. Since the cosmic density 
of the gravitino should be less than the dark matter density of the universe ^7j, we obtain 
the constraint 

Th < 3 X 10^ GeV ( "^'^ \ ' ( ^^"l for 7713/2 ~ 10"^ - 10 GeV, (17) 
V500GeV/ VlGeVy ^' ' v ; 

where we have omitted the logarithmic corrections. For 1 keV ^ 70.3/2 ^ 10^^ GeV, the 
upper-bound on Tr is of the order of 100 GeV, 

Tr < O(IOO) GeV for 7773/2 ^ 1 keV - 10"^ GeV. (18) 
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When the gravitino is hghter than 1 keV, no constraint comes from the cosmic density. 
However, such a hght gravitino behaves as warm or hot dark matter component and 
affects the power spectrum of the density fluctuations through free streaming. This may 
extend the bound (HBD to mg/s ~ 0(10) eV [22] • 

4 Gravitino production in inflaton decay and its cos- 
mological constraints 

In this section we first estimate the decay rate of the inflaton into a pair of the gravitinos 
and clarify a condition under which this decay channel becomes effective. Then we discuss 
cosmological constraints on such a decay for a broad range of the gravitino mass: ^3/2 = 
1 keV — 100 TeV, in order to show how severe this gravitino-oveproduction problem is. 

4.1 Inflaton decay into a pair of gravitinos 

The relevant interactions for the decay of an inflaton field into a pair of the gravitinos 
are |23 

o 

where ip^ is the gravitino field, and we have chosen the unitary gauge in the Einstein 

frame with the Planck units, Mp = 1. We have defined the total Kahler potential, G = 

K + ln where K and W are the Kahler potential and superpotential, respectively. 

Here and in what follows a subscript i denotes a derivative with respect to the field i, 

while a superscript is obtained by multiplying with g^^* , the inverse of the Kahler metric 

Qij* = Gij*. The SUSY breaking field z is such that it sets the cosmological constant to 

be zero, i.e., G^Gz — 3 and we assume that 2; is a singlet under any symmetries as 

in the gravity- mediated SUSY breaking models. In fact, the existence of a singlet (and 

^^Throughout this paper we assume that the _D-term potential is neghgible. In a broad class of the 
SUSY breaking models, the z field may not be the only field that has a sizable F-term, and \Gz\ could 
differ from \/3. However, the following arguments remain virtually intact as long as \Gz\ — if G^Gz 

decreases by one order of magnitude, the constraints on inflation models would become relaxed by the 
same amount. 
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elementary) field z with a nonzero F-term of 0{m3/2Mp) is a generic prediction of the 
gravity-mediated SUSY breaking models. This is true even in the case of the dynamical 
SUSY breaking [21], because the gauginos would become much lighter than squarks and 
sleptons otherwise j2S] Later we will give a comment on the case that z is charged 
under some symmetry as in the gauge-mediated and anomaly-mediated SUSY breaking 
models. 

It has been recently argued that the modulus and infiaton decays produce too much 
gravitinos through the above interaction j2Zll2ll2Hl- Taking account of the mixing between 
(f) and z, however, the effective coupling of the infiaton with the gravitinos is modified |29j . 
According to the detailed calculation of Ref. [201 , '^^ '^'^^Y have to replace with Q^^^ 
(the relation between these two is given in the next section). The real and imaginary 
components of the infiaton field have the same decay rate at the leading order [2Z1 ISDl '■ 

r,,,.r(^^2v-3/.)^f^^. (20) 

where we have assumed that the infiaton has a supersymmetric mass much larger than 
the gravitino mass: ^ ^3/2- Thus the decay rate is enhanced by the gravitino mass 
in the denominator, which comes from the longitudinal component of the gravitino as 
emphasized in Ref. |27j . 

It should be noted that the above expression for the decay rate cannot be applicable 
for H > 171^/2. The decay proceeds only if the Hubble parameter H is smaller than 
the gravitino mass, since the chirality flip of the gravitino forbids the decay to proceed 
otherwise. Intuitively, the gravitino is effectively massless as long as H > m^i2- 

We should clarify another important issue: what is the longitudinal component of the 
gravitino (i.e. goldstino) made of ? A similar issue was discussed in the context of the 
non-thermal 'gravitino' production during preheating [^, and it was concluded that the 
infiatino, instead of the gravitino in the low energy, is actually created [S3| Since the 

^•^Note, however, that it is possible, though comphcated, to generate a sizable gaugino mass by intro- 
ducing extra chiral superfields in the adjoint representation of the gauge group, rather than a singlet pBj. 

"^^^The decay can also be understood in terms of the goldstinos due to the equivalence theorem in 
supergravity 

"^^It should be noted, however, that the inflatinos produced during preheating may be partially con- 
verted to the gravitinos in the low energy, since is generically nonzero in the true minimum 2 (the 
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inflatino decays much earlier than the BBN epoch |31], the non-thermal 'gravitino' (ac- 
tually, infiatino) production turned out to be harmless. The reason is that the 'gravitino' 
production occurs in a rather early stage of the reheating just after the inflation ends, 
during which the energy stored in the inflationary sector significantly contributes to the 
total SUSY breaking. In our case, however, the situation is completely different; the decay 
into the gravitinos is effective, since we consider a cosmological epoch, H < m3/2, when 
the SUSY breaking contribution of the inflaton is subdominant. Thus, the gravitinos 
produced directly by the inflaton decay should coincide with those in the low energy. 

Let us now consider the implication of ()20j) . As we will see in the next section, 
the effective coupling Q^^^ is proportional to G,^, for such non-minimal interaction as 
(K/2)|0pz2; + h.c. in the Kahler potential. The auxiliary field Gcf, represents the fractional 
contribution of the inflaton to the SUSY breaking. One might suspect that (and 
therefore Q^^^^) should be zero in the vacuum and such a decay does not occur at all. 
However, as we will see in the next section, this is generically not true. To be sure, in 
many inflation models, the minimum of the inflaton potential preserves SUSY, as long 
as the inflaton sector is concerned. But, once we take account of the SUSY breaking 
sector, the minimum slightly shifts and non-vanishing is induced. This means that we 
need to consider the scalar potential including both the inflaton and the SUSY breaking 
sector field, in order to evaluate (and Q^^^^). Our next concern is how large Q^^"* 
can be. According to the general formula in single-inflaton models to be derived in 
Sec. El (see ()43|) ). it is at most ~ In fact, this is also true in the inflation 

models with multiple fields. Therefore the decay rate (j2(jp can be comparable to that 
obtained by the decay via Planck-suppressed dimension 5 operators. In other words, this 
direct gravitino production becomes important especially when the total decay rate of 
the inflaton is suppressed, i.e., the reheating temperature is low. Therefore the direct 
gravitino production via the interaction ()19|) is complementary to the thermal gravitino 
production which becomes more effective for higher Tr. Fig. ^ schematically shows this 
feature. This specific character enables us to put severe constraints on inflation models. 

inflation model adopted in Ref. \,VA\ has vanishing G^). This effect may further constrain the inflation 
models. 
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Y3/2 A 




direct^ ^ 



thermal 



Tr 



Figure 1: Dependence of the gravitino-to-entropy ratio on the reheating temperature T^. 
The sold line represents the abundance of the thermally produced gravitinos, while the 
dashed line corresponds to that directly produced by the infiaton decay. See © and 
in the text. 

4.2 Cosmological constraints on Q^^^^ 

In the following we assume that the reheating temperature satisfies the bounds from the 
thermally produced gravitinos discussed in Sec. 0*^. The reheating temperature Tr is 
related to the decay rate of the infiaton into the SM particles (and their superpartners) 



where counts the relativistic degrees of freedom and hereafter we set = 228.75. 
When the Hubble parameter becomes comparable to T^m, the infiaton decays. It is easy 
to see that H ~ Fsm ^ "^3/2 is reahzed at the decay, if the reheating temperature 
satisfies the bounds from the gravitinos produced by thermal scattering (i.e., (jlj) — (fTUjl . 
(fT^. and (Uni) — (fTSj)). Therefore the infiaton decay into the gravitinos is effective. The 

^^As pointed out in Ref. [20], the mixing between the inflaton and the SUSY breaking field may enhance 
the reheating temperature. Including the effect of the mixing, we take the reheating temperature as a 
free parameter throughout this paper. 

^''We have assumed r3/2 ^ Tsm, since the standard cosmology would be upset otherwise. 



by 




(21) 
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gravitino-to-entropy ratio is then given by *® 



3/2 — — -p, , ) 

where we have neglected the gravitino production from the thermal scattering. 

First let us consider the cosmological bound on the gravitino abundance for stable 
gravitinos of m3/2 ^ lOGeV. The gravitino abundance should not exceed the dark matter 
abundance; 

m3/2>^3/2 < ^^DM^ < 4.7 X 10-1° GeV, (23) 

where pc is the critical density, and we used ^dm^^^O.IS at 95% C.L. in the second 
inequality. Combining and (j^ . we obtain 

ier'l<3.2.:o-"(^)*(^)-(^)', m 

for Tr satisfying (fTTj) or (fTHj) . To further reduce this bound, we need to substitute the 
largest allowed value of Tr given by (fTTj) and (fTHj) . Then we arrive at 

< O,0-.«)(^)'(^)- ,25) 

for 7713/2 — 1 keV — 100 keV, and 

' ~ UoOGevj llGevj UoiOGevJ ^^^^ 

for m3/2 — 100 keV — 10 GeV. It should be noted that the constraints on Q^^^"* become 
severer for lower Tr, as clearly seen from (or Fig. Q). 

Next we consider unstable gravitinos. The gravitino abundance is severely constrained 
by BBN as discussed in Sec. 0] We can similarly derive the constraints on Q^^^^ from Q 
"i^^Here we assume that the entropy comes solely from the perturbative decay of the inflaton. 
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(jT^. and (1221): 
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- 30 TeV 


for 


"^3/2 
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- 1 TeV 


for 
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3 TeV 


for 


"23/2 




3 - 


10 TeV 



(27) 



for Bh — 10 ^. For 7/13/2 larger than 30(10) TeV, the constraint comes from the LSP 
abundance produced by the gravitino decay. Using Q, (fT^ . and (j^ . we obtain 



< 5 X 10" 



UOOGeV/ 



-1 



/ ^3/2 
UOOTeV; V 1010 GeV 



for 1713/2 — 30(10) — 100 TeV. In particular, this can be rewritten as 

-2 



< 2 X 10" 



1010 GeV 



(29) 



(30) 



for the anomaly-mediated SUSY breaking with the wino LSP, where we have used fllSj) . 

In Figs. 121 - in we show the upper bounds on Q^^^ together with predictions of new, 
hybrid, smooth hybrid, and chaotic inflation models to be derived in Sec. for repre- 
sentative values of the gravitino mass: 1713/2 = 1 GeV, 1 TeV, and 100 TeV, respectively. 
From these figures one can see that the bound is the severest in the case of ^3/2 = 1 TeV 
due to the strict BBN bounds. The bounds are slightly relaxed for either (much) heavier 
or lighter gravitino mass. Note that the constraints on inflation models do not change 
for the stable gravitinos with 1113/2 — 100 keV — 10 GeV, since both the upper bound and 
the actual value of Q^^^ in the vacuum are proportional to ^3/2 (cf. fj26|) and (021) )• The 
smooth hybrid inflation is excluded for a broad region of the gravitino mass, unless k (see 
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m3/2 = 1 GeV 




10° 
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new(multi) 
hybrid 
smooth hyb. 
chaotic 



108 10IO 10I2 1014 10I6 



[GeV] 



Figure 2: Upper bound on the (effective) auxiliary field of the inflaton as a function 
of the inflaton mass m^, with 7713/2 = 1 GeV. We set = 500 GeV. is set to be the 
largest allowed value, and the bound becomes severer for lower Tr. The typical values of 
Q^^^ and for the multi-field new, hybrid, smooth hybrid, and chaotic inflation models 
with K = 1 are also shown. The chaotic inflation can avoid this bound by assuming Z2 
symmetry (see the text for details). 

the next section for the definition) is suppressed. Similarly, for n ~ 0(1), a significant 
fraction of the parameter space in the hybrid inflation model is excluded, and in partic- 
ular, it is almost excluded for 7773/2 = ITeV, while the new inflation is on the verge of. 
Even though the constraints on the hybrid inflation model seems to be relaxed for smaller 
777,^, it is then somewhat disfavored by WMAP three year data [20j since the predicted 
spectral index approaches to unity. The chaotic inflation model is also excluded unless k 
is suppressed due to some symmetry. 



5 Explicit calculation of for several inflation mod- 



In estimating the effective coupling of the inflaton (j) with the gravitino, the mixings with 
the SUSY breaking sector field is important, as pointed in Ref. [21] for specific cases. In 
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X : chaotic 
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Figure 3: Same as Fig. |21 except for 7/13/2 = 1 TeV. The typical values of and 
for the single-field new inflation model with k, = 1 are also plotted. The solid and dashed 
lines are for the hadronic branching ratio = 1 and 10"'^, respectively. 

fact one can rigorously estimate the coupling in a rather generic way ^30| . In this section, 
we would like to show how to obtain based on the argument of Ref. |3Uj . 

The point is that the inflaton field does not coincide with the mass eigenstate after 
inflation due to the mixings with the SUSY breaking sector field z. There are three sources 
for the mixings: (i) kinetic terms (or equivalently, the Kahler metric); (ii) non-analytic 
(NA) mass terms; (iii) analytic (A) mass terms. Although the mixing in the Kahler metric 
can be important for the reheating processes jHOj, we neglect it here since it does not affect 
the coupling with the gravitinos. In the following we focus on the mixings in the mass 
terms. 

5.1 Single- field inflation model 

Let us first consider a single- field inflation model, with the Kaler metric Qij* = 6ij. In the 
Einstein frame, the SUGRA Lagrangian contains the scalar potential, V = e'^iG'^Gi — 3). 
The non-analytic (NA) and analytic (A) mass terms are written as 




{v iGk^ j^G'^ — Rij*ki*G^G^* + Qij*^ , 



(31) 
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m3/2 = 100 TeV 




+ : new(single) 

: new(multi) 

: hybrid 

: smooth hyb. 

X : chaotic 



[GeV] 

Figure 4: Same as Fig. |21 except for 7/13/2 = 100 TeV. 



M^. = Mf, = = (V.G, + V,G. + G'V.V.G,) , (32) 

respectively, where we have assumed the vanishing cosmo logical constant, C'Gi = 3, and 
used the potential minimization condition, G'^SI kGi+Gk = in the vacuum. The gravitino 
mass is given by 7713/2 = (e'^^'^^. Here Rij*ke* is the curvature of the Kahler manifold, 
defined by Rij*ki* = gij*ki' — 9^^* 9mj'i*gn*ik- Also the covariant derivative of Gi is defined 
by ViGj = Gij — T^jGk, where the connection, T^j = g^^* gije*, and Vkgij* = is satisfied. 

We assume that the inflaton is heavy due to a large supersymmetric mass, = 
|gG/2y^^^| ^ 7713/2, and that M^^ dominates over the other elements of the mass terms. 
Then the NA mass terms can be diagonalized by the following transformation: 

$ = + ez, 

Z = z-e*(j), (33) 

where e represents the mixing angle. Here we have assumed |e| <^ 1 and neglected those 
terms of O(e^). Since M?t dominates over the other components in the mass matrix, the 
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mixing angle is given by the ratio of M^^ to the off-diagonal component: 

e - T#- (34) 

As emphasized in Ref. the NA mass eigenstates {^,Z) do not necessarily coincide 
with the true mass eigenstates. In fact, the analytic mass terms generically provide further 
mixing between and z^. The true mass eigenstates are therefore 

$ = (j) + ez + ez\ (35) 
Z = z-e*(f)-i(j)\ (36) 

where the mixing angle e, which is assumed to be much smaller than unity, is given by 

Mt 

(37) 



ML 

<p<p 

Below we show that the coupling with the gravitinos is suppressed in the NA mass eigen- 
states, but it is not the case in the true mass eigenstates if the Kahler potential is non- 
minimal. 

In the NA mass eigenstates, the off-diagonal element of the non-analytic mass term is 
zero by definition: 

Mlz = (V^G^V^G^ + V^G^V^G^ - i?$^,,*G'G^*) = 0, (38) 
which leads to 

V,G, ^ (39) 

where we have used |V$G$| ^ \^zGz\- On the other hand, the potential minimization 
condition for $ reads 

GiV$G$ + GzV^Gz + G$ = 0, (40) 

which can be solved for G$: 

G, ^ -^Gz. (41) 
Substituting into (jH)), we arrive at 

|G,| ^ 3v^^*f4, (42) 

I V(j)(jr<j, I 
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where we have used \Gz\ = \G^\ — v^. Thus G$ is always proportional to 'rn\i2/m?^ <C 1. 
For the minimal Kahler potential, Gq, is exactly zero in this basis. In the mass-eigenstate 
basis, therefore, the effective coupling of the infiaton $ with the gravitinos dominantly 
comes from the mixing in the analytic mass terms: 

\Gf \ ^ \eG^^?,\g-,j'^. (43) 

For instance, let us consider 5K = {l/2)K,\(f)\'^zz + h.c, which is expected to be present 
if z is singlet under any symmetries. For the non-minimal Kahler potential, the effective 
coupling becomes 

l^ri - 3.(0)^, (44) 

where {(f)) denotes the vacuum expectation value (VEV) of (f). Therefore, Q^^^ is propor- 
tional to G^ ~ (0) Tn^i2/Tn^ 2j, for this choice of the interaction between cf) and z. 

Here let us comment on the case that the mass of the SUSY breaking field Z, mz-, is 
larger than due to the non-SUSY mass term. Such situation may be realized in the 
dynamical SUSY breaking models [21]. Then the l^^^^'*] becomes of O((0) m^/2/^(i)) even 
if the Kahler potential is minimal [211 IHO] • To be conservative, however, we assume that 

> irtz in the following discussion. 

As a concrete example, here we study the new inflation model jHSlEElEZI- In the new 
inflation model, the Kahler potential and superpotential of the infiaton sector are written 
as 

K(0,0t) = 101^+^-'^ 



W{<P) = v^<P ^0"+^ (45) 



4' 

n + 1 

where the observed density fluctuations are explained for t> = 4 x 10^^(0.1/(7)^/^ and 
fc^ 0.03 in the case of n = 4 [37]. After inflation, the infiaton takes the expectation 
value (0) ~ {v^ jg)^!"^. In this model the infiaton mass is given by ~ nv^ j (0), and 
the gravitino mass is related to v as ^3/2 — nv^ (0) /(n + 1), since the infiaton induces 



18 



the spontaneous breaking of the i?-symmetry. Thus, ()4H|1 leads to *^ 



For the interaction 5K = (l/2)/t|0p2;z + h.c, this becomes 



eri-r^tejf^V. (46) 



3_ 

t2 \ n 



n + 1 \ g 



In the case of n = 4, 1^^^''^'' | ^ 8 x 10 and ~ 4 x 10^ GeV for 7713/2 = 1 TeV, while 
\g^^^^\ ~ 8 X lO^^K and ~ 2 x 10^° GeV for mg/a = 100 TeV. Note that mg/a < 1 TeV 
cannot be realized unless g ^ 1. We plot these results with k = 1 in Figs.|niand|31 We can 
see that the new inflation model is on the verge of being excluded for 7713/2 = 1 TeV ^^'^j 
while it is close to but slightly below the bound for 7/13/2 = 100 TeV. This single- field new 
inflation model will be discussed in detail in Ref. 

5.2 Multiple-field infiation model 

Next we consider an infiation model with multiple fields, for which the formula ()4H|1 cannot 
be simply applied as it is. Although we generically need to evaluate Q^^^^ for each infiation 
model, there is an important class of models described by the following superpotential: 

= (48) 

where f{ip) is a function of ip. The potential minimum in the global SUSY limit is located 
at 

(0) = 0, 

W = ^0, (49) 

where ipo satisfies f{ipo) = 0. Note that the true minimum is slightly displaced from ()49|) . 

once the SUSY breaking field is taken into account [21 HO] . 

'^^The relation 1)43(1 remains virtually unchanged in the presence of the quartic coupling in the Kahler 
potential. 

'^^"it may survives if k ~ 10^^ as suggested in the large-cutoff SUGRA [^. We thank M. Ibe and Y. 
Shinbara for useful discussion. 
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For instance, the above class of the models includes a new inflation model [H] and a 
hybrid inflation model |13 113 HI] , described by 



W{<P^^) = (50) 



where determines the inflation energy scale and M is an effective cut-off scale. In the 
new inflation model %p plays a role of the inflaton, while is the inflaton in the hybrid 
inflation model. 

The inflaton fields and have almost the same masses, 



(51) 



which are assumed to be much larger than the gravitino mass. It should be noted that 
and ip (and/or ip'^) almost maximally mix with each other to form the mass eigenstates 
due to the almost degenerate masses. To see this let us take the NA mass-eigenstate basis 
($, Z) in which the non-analytic mass matrix is diagonalized except for 4> — ip mixing. 
The difference between the diagonal components of the non-analytic mass matrix is small: 
|M||, — M^^l = 0(7713^2)) while the off-diagonal component in the analytic mass matrix is 
relatively large: = 0(m3/2^(/,) resulting in the almost maximal mixing between 

and ipK This mixing is effective at the inflaton decay, since the Hubble parameter at the 
decay should be (much) smaller than O {1713/2) to satisfy the bounds from the thermally 
produced gravitinos. However, since the mixing is due to the specific character of (j48j) 
and it occurs within the inflaton sector, we leave it for a moment. Then we can similarly 
show that the auxiliary fields G$ and Gq, are proportional to ^3/2/'^0 fhe NA mass 
eigenstates (<I',\E',Z). Therefore the effective couplings with the gravitinos arise mainly 
from the mixings in the analytic mass terms, as in the single- field inflation *^^: 

l^^^^^l - S\g^J^. (52) 



'^^^In addition, the off-diagonal component in the non- analytic mass matrix as well can be as large as 
-^^$\i> = 0{m^/2ni(j,) if f'i'ipo) ^ /'('0o)/^o, and the mixing is almost maximal in this case too. 

#^^Note that the dependence of the right-handed side on (f) and ip originates from the SUSY mass (|51|l . 
which is peculiar to the form of the superpotential jSUJ. 
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For such interactions as 5K = {n/2)\%p\^zz + (/t/2)|0p2;z + h.c, we have 

l^ri ^ 3^(0)^. (53) 

Therefore l^^'^^^l is suppressed compared to l^^'^^l if ^ (0) as in the case of ^E^. 

The true mass eigenstates are obtained after taking account of the (almost) maximal 
mixing between and '0('0^) discussed above: 

V± - ^ , (54) 



where we have omitted the relatively small mixings with z for simplicity, but they are 
included in the definition of ip±. Fo 
the gravitinos are roughly given by 



included in the definition of ip±. For |^$''^''| ^ |^*^^|, the effective couplings of (p± with 



l^if I - -^\Qf\ (55) 



5.2.1 New inflation model 



The new inflation discussed in Sec. 15. II is also realized for 

K = \<P\' + + 'j\<P\' + kM'm' + ji^i', 

W = (Piv^-gij^), (56) 

in which the infiaton is ip, while stays at the origin during and after inflation If 
one defines k = k2 — 1, the scalar potential for the infiaton ip becomes the same as the 
single-field new inflation model, although the gravitino mass is not related to the infiaton 
parameters. After the inflation ends, the energy of the universe is dominated by the 
oscillation energy of tjj Although q!^^^ is suppressed compared to G^^\ the effective 
coupling to the gravitinos is given by (j55|) . since and ip almost maximally mixes with 

i^^^If one introduces a constant term in the superpotential, the (p shifts from the origin. 

#i4rpj^g tachyonic preheating |45l I4f)j is known to occur in this model, and if it occurs, the homogeneous 
mode of the infiaton ip disappears soon and the excited tp particles are produced. This instability itself 
does not relax the gravitino-overproduction problem, since these iJj particles will decay perturbatively 
into the SM particles and their superpartners. Further, if the 4' particles are relativistic, the decay is 
delayed, making the problem even worse. 
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each other in the vacuum. Thus the constraint on this model is comparable to that on 
the single-field new infiation. For the non-minimal coupling 5K = {l/2)n\ip\'^zz + h.c, 
the effective coupling to the gravitinos is given by 

We plot the value of for g = IQ-^ - 1 and k = 10"^ - lO"^-^ with the e-folding 

number = 50 in Figs. |2l El and 0] Thus the (multi-field) new infiation model is on the 
verge of being excluded, if k is order unity. 

5.2.2 Hybrid infiation model 

The hybrid infiation model contains two kinds of superfields: one is which plays a role 
of infiaton and the others are waterfall fields and |42| |331 Hlj- After infiation ends, 
(f) as well as iplip) oscillates around the potential minimum and dominates the universe 
until the reheating. 

The superpotential W{(j),ip,4') for the infiaton sector is 

W{<P,i;,^) = <Pifx'-\i,^), (58) 

where ip ^^^^ assumed to be charged under f/(l) gauge symmetry. Here A is a 

coupling constant and fi is the infiation energy scale. The potential minimum is located 
at (0) = and {^jj) = (?/') = /i/-\/A in the SUSY limit. For a successful infiation, /i and A 
are related as ~ 2 x lO^^A^/^ for A ^ 10"^, and /i ~ 2 x lO^^A^/® for A ^ 10^^. Moreover, 
in this type of hybrid infiation there exists a problem of cosmic string formation because 
ip and ip have U{1) gauge charges. To avoid the problem the coupling A should be small 
as, A ~ 10-4 [47]. 

Due to the D-term potential one linear combination of ip and ip, given by ip^~^ = 
{ip — ^)/\/2, has a large mass of ~ (7 {ip) {g denotes the gauge coupling), while the other, 
ip^^^ = (ip + ■?/;) /v^ has a mass equal to that of (p: m^(+) = = ^/2\{ip). It is the latter 
that (almost) maximally mixes with (p to form mass eigenstates. Since the form of the 
superpotential is almost identical to (j^ . it is straightforward to extend the results 
and to obtain 
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for the non-minimal coupling 5K = (l/2)K,\'ip^'^^\'^zz + h.c. Note that VEV of ijj^^^ is 
equal to \/2{ip). 

For A ~ 10-1 - 10-^ [IH: we obtain ~ 8 x 10"^ - 1 x 10-^ |^Jf)| ~ 9 x 
10-i5K(m3/2/l TeV) - 9 x 10-^^K{ms/2/l TeV) and ~ 8 x 10^^ - 1 x 10^° GeV. From 
Fig. El one can see the hybrid inflation model is almost excluded by the gravitino over- 
production for 1713/2 = 1 TeV, if k is order unity. For m3/2 = 1 GeV and 100 TeV, the 
constraints become slightly mild, but a significant fraction of the parameter space is 
still excluded (see Figs. El and E)). Although the constraints on become relaxed 

for smaller rrifp (i.e., smaller A^10~^), it is then somewhat disfavored by the WMAP 
data [20j since the density fluctuation becomes almost scale-invariant. 

Next let us consider a smooth hybrid inflation model |49^, which predicts the scalar 
spectral index as Ug — 0.97, which is slightly smaller than the simple hybrid inflation 
model. The superpotential of the inflaton sector is 

WM,^)=4>{f^'-^)- (60) 

The VEVs of tfj and tp are given by {ip) = (^ip^ = {jj.M'^'^Y^^ ^ and we assume that ip = ip 
always holds due to the additional D-term potential. Then one of the combination, 
mixes with 0, and the effective coupling with the gravitinos is given 

by 

I A,(V2(^))!^, (61) 

for 5K = (1 /2) k\^^~^'' zz + h.c. Here we have defined m^j, = \/2nn'^/ {ip). The constraint 
on the model is more or less similar to that on the hybrid inflation model. In fact, for 
= 2 we obtain /i ~ 4 x 10"^ - 9 x 10-^ ~ 2 x 10-i3fi:(m3/2/l TeV) - 4 x 

10-i'^K(m3/2/l TeV) and ~ 1 x 10^^ - 6 x 10^^ GeV. From Figs. El- Hone can see that 
the smooth hybrid inflation model is excluded for a broad range of m3/2 for k = 0(1). 

Lastly let us comment on the D-term inflation model [^0] , in which one of the waterfall 
fields, ip-, obtains a large VEV of ~ 10^^ GeV. The ip- field can decay into a pair of the 
gravitino if there is a coupling like \ip-\'^zz/2 + h..c. in the Kahler potential. However, after 
inflation, the universe is dominated by the two fields: one is the ip- field and the other is 



23 



the inflaton, 5*. The resultant gravitino abundance thus depends on both the reheating 
processes of these two fields and the relative portion of the energy in each field |51j . 
Therefore we cannot put a rigorous bound on the D-term infiation model. 

5.2.3 Chaotic inflation model 

A chaotic infiation is realized in SUGRA, based on a Nambu-Goldstone-like shift 
symmetry of the infiaton chiral multiplet El] . Namely, we assume that the Kahler 
potential K{(f), 0^) is invariant under the shift of 0, 



+ iA, 



(62) 



where A is a dimensionless real parameter. Thus, the Kahler potential is a function of 
+ 0^; K{(f), 0"!') = K{(j) + 0"^) = c (0 + 0"'") + 1(0 + 0^)^ + ■ ■ where c is a real constant 
and must be smaller than 0(1) for a successful infiation. We will identify its imaginary 
part with the infiaton field if = \^lm[(f)]. Moreover, we introduce a small breaking term 
of the shift symmetry in the superpotential in order for the infiaton ip to have a potential: 



W{(f), ip) = m(f)ip, 



(63) 



where we introduced a new chiral multiplet ip, and m ~ lO^^GeV determines the infiaton 
mass. 

The scalar potential is given by 



2„K 



V{r], ip,4>) = m e 



1^1' l + 2h + 



V2, 



7] + {V + 



+ m 



V2, 



with 



(64) 



(65) 



where we have assumed the minimal Kahler potential for ip, and defined r] = A/2Re[0]. 
Note that rj and cannot be larger than the Planck scale, due to the prefactor . On 
the other hand, ip can be larger than the Planck scale [S3] , since cp does not appear in K. 
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For if ^ 1, r] acquires the mass comparable to the Hubble parameter and quickly settles 
down to the minimum, 77 ~ —c/ -\/2- Then the scalar potential during inflation is given by 

V{r],(f,ip) ~ ^m'^ip'^ + ■m'^\ip\'^. (66) 

For ^ 1 and < 1, the ip fleld dominates the potential and the chaotic inflation takes 
place (for details see Refs [5^1 EH)- 

The effective auxiliary fleld of ip is given by 

m Sg^,, — = 3k — '-, (67) 
m m 

where we have assumed the non-minimal coupling SK = (1/2)k(</) + (f)^)zz + h.c. in the 
second equality. This Kahler potential is invariant under the shift symmetry (jU^ . Note 
that l^^'^^^l is suppressed for e.g., 6K = {l/2)k\il)\^zz + h.c. due to {i})) <^ 1. Taking 
account of the mixing between and tp^ the effective coupling with the gravitinos is 
given by 

(68) 

It is worth noting that both real and imaginary components of </> can decay into a pair 
of the gravitinos via the mixings with z and if). One might suspect that it is only the 
real component of (j) that can decay into the gravitinos, since the shift symmetry dictates 
that the only real component (0 + 0^) appears in the Kahler potential. However, it is not 
surprising that this is not the case, since the enhanced decay amplitude is proportional 
to powers of the large SUSY mass m that explicitly violates the shift symmetry. 

We plot the result (|UHj) with k = 1 in Figs. El El and El Although the coupling is too 
large if k = 0(1), it should be noted that in this chaotic inflation model we can realize 
~ by assuming an approximate Z2 symmetry. Therefore, the new gravitino 
problem does not exist in this case. A detailed discussion on the chaotic inflation model 
will be given in [HH] . 

6 Conclusions 

Throughout this paper we have assumed no entropy production late after the reheating of 
inflation. We briefly discuss potential problems when a late-time entropy production 
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occurs. First of all, the cosmological constraints on the reheating temperature shown in 
Sec. ini would be relaxed and the Hubble parameter at the inflaton-decay time does not 
necessarily satisfy the condition H < 777,3/2 for the formula (j^^ to be applicable. Thus, the 
cosmological constraints on Q^^^ would become milder On the other hand, we must 
be careful about the gravitino production in decay processes of the field X responsible for 
the late-time reheating. One may have a similar stringent constraint on Qx^^ . An obvious 
way to induce late-time entropy production avoiding the problem is to assume the late- 
time decay of a scalar field with a mass smaller than 27773/2. In addition, there is another 
interesting example that is free from the problem. Consider that the scalar partner of a 
right-handed neutrino possesses a large value during the inflation. If the value is at the 
Planck scale and its decay rate is small, the scalar dominates the universe before its 
decay. Thus, the decay of the scalar A^ can produce entropy and dilute the abundance of 
the relic gravitino. The crucial point here is that the scalar A^ does not decay into a pair 
of gravitinos due to the matter (or lepton-number) parity conservation. In other word, 
A^ does not mix with the SUSY breaking field. Thus, this decay process is free from the 
gravitino-overproduction problem. Furthermore, the decay of the scalar A^ may generate 
the baryon asymmetry of the universe p7j through the leptogenesis |[5H] . 

Another even manifest solution to the gravitino overproduction problem is to assume 
the gravitino mass 7773/2 < 0(10) eV [221- In this case, the produced gravitinos get into 
thermal equilibrium due to relatively strong interactions with the standard-model parti- 
cles, and such light gravitinos are cosmologically harmless. 

Let us comment on another decay mode induced by non-minimal couplings between 

the infiaton and the SUSY breaking field z. From ()43|) and ()52|). the gravitino production 

rate is proportional to \g^zz\'^- If 94>zz is nonzero, the infiaton (p can also decay into the 

SUSY breaking field z jSO], and the partial decay rate is comparable to that into the 

gravitinos. As noted in Ref. thus produced z may cause a cosmological problem at 

most as severe as that induced by the gravitinos. Therefore including the effect of the z 

T^^^Even if the reheating temperature is higher than the cosmological bounds discussed in Sec. 15.11 
the direct gravitino production by infiaton decays can occur and the formula (|20() is applicable as 
long as the condition H < 7713/2 is satisfied at the decay time. In this case, we must consider the 

direct production with a great caution, since it dominates over the thermal production if > 
10-^3(Tij/10iOGeV)(10i3GeV/m^)2(m3/2/lTeV). 
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production may make the problem only a few times worse, and our discussion remains 
qualitatively unchanged. 

In this paper we have shown that an inflation model generically leads to the grav- 
itino overproduction, which can jeopardize the successful standard cosmology. We have 
explicitly calculated the effective auxiliary fleld Q^^\ which is an important parameter 
to determine the gravitino abundance, for several inflation models. The new inflation is 
on the verge of being excluded, while the (smooth) hybrid inflation model is excluded if 
K = 0{1). To put it differently, the coefficient of the non- minimal coupling in the Kahler 
potential, k, must be suppressed especially in (smooth) the hybrid inflation model. We 
show the constraints on k for the inflation model we studied so far in Figs. El-IZI As long 
as the SUSY breaking fleld z is singlet, there is no reason that k should be suppressed. 
Therefore those inflation models required to have k ^ 1 involve severe flne-tunings on 
the non-renormalizable interactions with the SUSY breaking fleld, which makes either the 
inflation models or the SUSY breaking models containing the singlet z (with Gz = 0(1)) 
strongly disfavored. We stress again that the existence of such a singlet fleld is required in 
the gravity-mediated SUSY breaking, in order to give the SM gauginos a mass comparable 
to the squark and slepton masses. One of the most attractive ways to get around this new 
gravitino problem is to postulate a symmetry of the inflaton, which is preserved at the 
vacuum, to forbid the mixing with the SUSY breaking fleld. Among the known models, 
such a chaotic inflation model can avoid the potential gravitino overproduction problem 
by assuming Z2 symmetry. Another is to assign some symmetry on the SUSY breaking 
fleld z as in the gauge-mediated and anomaly-mediated 3j SUSY breaking models. So 
far we have assumed that z is singlet under any symmetries as in the gravity-mediated 
SUSY breaking models. If the SUSY breaking fleld z is not a singlet, and the non- minimal 
coupling like 6K = K/2\(f)\'^zz + h.c. can be suppressed. It should be noted however that 
the mixing between and z may induce other cosmological problems jHO] even if z is 
charged under some symmetry and/or its VEV is suppressed. 

Although we have briefly discussed various (typical) inflation models, it should be 
stressed that the gravitino-overproduction problem is common to all the inflation models 
in SUGRA. Thus, in inflation model building, one must always check whether an inflation 
model under consideration satisfles the bound. 
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Figure 5: Upperbounds on k as a function of the inflaton mass m^, with 1713/2 = 1 GeV, 
for the multi-field new, hybrid, and smooth hybrid, and chaotic inflation models are also 
shown. See the text for the definition of k in each model. 
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Figure 6: Same as Fig. El except for = 1 TeV. The bound on the single-field new 
infiation model is also plotted. 
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Figure 7: Same as Fig. IHl except for 1713/2 = 100 TeV. 
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